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Abstract 

We report results from a high resolution numerical study of fluid particles transported by 
a fully developed turbulent flow at Rx = 280. Single particle trajectories were followed for 
a time range spanning more than three decades, from less than a tenth of the Kolmogorov 
time-scale up to one large-eddy turnover time. We present results concerning acceleration 
statistics and the statistics of trapping by vortex filaments conditioned to the local values of 
vorticity and enstrophy. We distinguish two different behaviors between the joint statistics of 
vorticity and centripetal acceleration or vorticity and longitudinal acceleration. 

Lagrangian statistic of particles advected by a turbulent velocity field, u(x,t), is important 
both for theoretical implications [1] and for applications, such as the development of phenomeno- 
logical and stochastic models for turbulent mixing [2,3]. Recently, important advances in experi- 
mental techniques for measuring Lagrangian turbulent statistics [4-8] have been achieved. Direct 
numerical simulations (DNS) also offer very high accuracy albeit at a slightly lower Reynolds 
number [9-12]. We analyze Lagrangian data obtained from a recent DNS of forced homogeneous 
isotropic turbulence [17,18], performed on 512 3 and 1024 3 cubic lattices with Reynolds numbers 
up to i?A ~ 280. The Navier-Stokes equations were integrated using fully de-aliased pseudo- 
spectral methods for a total time T ~ Tl. Two millions of Lagrangian particles (passive tracers) 
were injected into the flow once a statistically stationary velocity field had been obtained. The 
positions and velocities of the particles were stored at a sampling rate of 0.07t v . The velocity of 
the Lagrangian particles was obtained using linear interpolation of the Eulerian field. Acceleration 
was measured both as the derivative of the particle velocity and by direct computation from all 
three forces acting on the particle (i.e. pressure gradients, viscous forces and large scale forc- 
ing): the two measurements were found to be in very good agreement (with the latter being less 
noisy) . Finally, the flow was forced by keeping the total energy constant in each of the first two 
wavenumber shells. For more details on the simulation, see [17-19]. Recently much attention 
focused on the statistics of acceleration with and without conditioning on the local structure of 
the flow [4, 8, 16, 18, 19]. Some phenomenological description of such statistics using multifractal 
[14, 18] or quasi-equilibrium distribution have been also proposed [13] (for a critical summary of 
these attempts see [15]). In this paper we concentrate mainly on the statistics of trapping events, 
i.e. those cases when the particle is captured inside a vortex filaments for a time lag consider- 
ably larger than the Kolmogorov eddy turn over time, t v [4, 17]. These events contribute to the 
statistics of the particle acceleration, V(a), with extremely intense values, up to 80 times the root 
mean squares acceleration, at the Reynolds number here investigated (see Fig. 1). In previous 
analysis [17-19] we have shown that the trapping in vortex filaments is responsible for a strong 
deterioration of scaling properties of Lagrangian Structure functions for time lags of the order 
of t v . In a previous publication [19] we shown that trapping into vortex filaments leads to very 
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different dynamical behavior for centripetal and longitudinal acceleration (see also [16]). Being the 
latter highly oscillating in time, while the former almost constant and with very high amplitude 
(see Fig. 1). The big difference between the dynamical properties of longitudinal and centripetal 
acceleration is already detected by inspecting their temporal correlations: 



C c (r) = (at(t)a x c (t + r)) d(r) = (a? (t)af (t + r)> 



(1) 



where we focus only on one component of the acceleration correlation for simplicity. In Fig. 2 
we present the DNS results for these two quantities. As one can see, while the longitudinal 
acceleration decreases with a characteristic time which is comparable with the Kolmogorov time 
r n , the centripetal one has a much slower time decay. Indeed, as shown in the inset, the centripetal 
correlation, C c (t), possess two different decays [4,5,16]. The first one, for short time lags (up 
to 4 -7- <3Tji) is comparable with the one of the longitudinal correlation. The second decay, with 
a characteristic time of the order of 10 -j- 15r„ is established for longer time lags. This second 
exponential tail is the signature of the strong persistence inside vortex filaments. In order to better 
quantify this trapping events, we present in this paper some results on the statistical properties 
of the longitudinal and centripetal acceleration conditioned to the local properties of the vorticity 
field, u>(x) and of the enstrophy, f2(aj). The presence of persistent intense vorticity structure at 
dissipative scales may have a strong feedback on the energy cascade mechanism. It is not clear if 
and how these structures may also influence the inertial range physics through non-local effects 
(filaments are quite elongated in one direction). Other studies similar to the one presented here, 
but focused on the joint probability distribution of the energy transfer and the vorticity field may 
help in clarifying this important issue. 

1 Joint statistics 

We want to study the correlation between the probability to observe a large value in the accel- 
eration fluctuations and the local Eulerian structure of the flow. To do that we first define the 
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Figure 1: (Left panel). Log- linear plot of the acceleration pdf, V (a). The crosses are the DNS data, 
the solid black line is the multifractal prediction [18] and the green line is the K41 prediction. The 
statistical uncertainty in the pdf is quantified by assuming that fluctuations grow proportional to 
the square root of the number of events. Inset: we show a V(a) to check the statistical convergence 
up to fourth order moments, the continuous line is the multifractal prediction presented in [18]. 
(Right panel). We show, in natural units, the behavior of one component of the centripetal and 
of the longitudinal acceleration. Notice the strong sign persistence of the centripetal acceleration 
with respect to the longitudinal one. 
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Figure 2: Correlation functions of the centripetal (+, red curve) and longitudinal (*, green curve) 
accelerations. Notice the much slower decay showed by the centripetal component. Inset: Log- 
linear plot of the centripetal correlation function, C c (r), with superposed the two exponential 
behavior for short time lags, exp(— r/(l.lr^)) and for large time lags, exp(— r/(15r ?) )), obtained 
with a best fit at short and long times respectively. 



longitudinal and centripetal instantaneous acceleration to be 

a; = (a ■ v)i> 

and 

a c = a x v 

respectively; where v(t) is the particle velocity unit vector, v(t) — u(r(t),t), and u(x,t) is the 
Eulerian velocity field. Similarly, we are interested to the square of the antisymmetric part of the 
stress tensor, the enstrophy, Sl(x): 

n{x)= l -Y,^- d ^) 2 =^ ( 2 ) 

y 

where ui is the vorticity. In a purely circular motion of radius r, around a vortex filament with 
vorticity, u>, we expect that the centripetal acceleration can be expressed as: 

|o c | oc u) 2 r - Qr (3) 

which is a direct links between the local properties of the enstrophy and the statistics of the 
centripetal acceleration. The above argument does not fix the 0(1) proportionality prefactor, 
which depends on the exact shape of the vortex filament. In order to probe how much the rare, but 
intense, events characterizing the tails of the acceleration statistics are indeed caused by trapping 
in regions of high vorticity we studied the joint probability densities of centripetal acceleration and 



3 



enstrophy, compared with the joint probability density of longitudinal acceleration and enstrophy 
(for recent similar investigation in experimental and numerical data, see [16]): 

7>(log(a c ),log(fi)), 7>(log(a z ),log(fi)) (4) 

where we concentrate directly with the logarithm of the variable to better appreciate extreme 
events. Let us first show in Fig. 3 the probability density function of the three -unconditioned- 
quantities: 

P(log(a c )), P(log(fi))), VQogfa)). 

According to relation (3) one would expect that if intense tails in the centripetal acceleration are 
due to the vortex trapping at Kolmogorov scale r], then the two PDF, "P(log(|a c |)) and V(log(Cl)) 
should have the same behavior for large values after the rescaling |a c |/?7 is made. This is what 
we showed in Fig. 3 where one can see that the two tails corresponding to V{\og(\a c \/rj)) and 
to 7 3 (log(il)) almost perfectly superpose ( within 0(1) prefactors which are out of control with 
dimensional arguments). 
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Figure 3: Logarithm of the Probability density function of the instantaneous centripetal accel- 
eration, V(\og(\a c \/ri)) (+, red curve), longitudinal acceleration V(log(\ai\/ri)) (*, blue curve) 
and of the enstrophy, V(log(Q)) (x, green curve). Acceleration amplitudes are rescaled with the 
Kolmogorov scale r\ to test relation (3). 

However, from the unconditional PDF shown in Fig. 3 there are not strong signs distinguishing 
the dynamical properties of centripetal and longitudinal acceleration. Similarly, the joint proba- 
bility densities (4) plotted in Fig. 4 do not show any quantitative differences between the statistics 
of the centripetal and longitudinal acceleration. This is quite natural, because the motion of a 
particle in a turbulent field will be characterized by different accelerations and decelerations, not 
necessarily associated with spiraling motion (on average the mean value of the acceleration will 
be zero) . The distinguishing character of vortex trapping is the strong persistence of the direction 
-and amplitude- of the centripetal acceleration for time lags much larger than t v , at variance with 
what happens to the longitudinal acceleration, see right panel of Fig. 1. 

To make this statement quantitative, we have studied the running average of the centripetal 
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Figure 4: Joint probability distribution function for one component of the centripetal acceleration 
and enstrophy, V(log(a^), log(fi)) (Left panel), and of longitudinal acceleration and enstrophy 
VQog(af ), log(fi)) (Right panel). Notice that the two shapes are almost indistinguishable. 

and longitudinal acceleration, over a time window varying up to IOt^ [19]: 

-, /.t+A/2 

a?{t) = - / dt'a c (t'); (5) 

A Jt-A/2 



We expect that the pdfs of the averaged centripetal and longitudinal acceleration will behave very 
differently at increasing the window size, A. In particular, the strong persistence of the centripetal 
acceleration up to IOtv, suggests that the joint PDF of centripetal acceleration and enstrophy, 
P(]og(af),log(n)) should remain almost unchanged at varying A, while the longitudinal one 
V(\og(a^), log(f2)) should experience a strong depiction, of events with simultaneous intense values 
of acceleration and vorticity. 

A way to offer a quantitative measurement of the correlation between the two variables entering 
the joint PDF, V{x,y), is to plot the residual value, 1Z(x,y) = V{x,y) — V(x)V(y), obtained by 
subtracting the probability density in the case of fully uncorrelated variables, V(x)V{y). This is 
what we show in Figs. (5-6) where we compare the residual value for the averaged centripetal 
acceleration and the local enstrophy: 




(6) 



ft(log(a*),log(n)) =P(log(a^),log(n))-P(log(a^))P(log(0)) 



(7) 



and the same quantity for the longitudinal one: 



ft(log(af ),log(ft)) = P(log(af),log(0)) - 7>(log(af ))P(log(fi)) 



(8) 
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Figure 5: (Left panel). Residual value, lZ(log(a^) , \og(Q)) for the centripetal acceleration and the 
local enstrophy without any average in time (i.e. A = 0). Acceleration is plotted on the y axis 
while the enstrophy is plotted on the x axis. (Right panel). The same but for the longitudinal 
acceleration, 7£(log(a^), log(fi)). The isoline is drawn for a fixed value of 1Z, to guide the eyes. 
Notice the strong correlation between intense vorticity events and intense centripetal and longitu- 
dinal acceleration (bottom right corner). Without averaging in time, no clear distinction between 
longitudinal and centripetal signal can be detected. As expected, the distinction arises only when 
persistent effects are investigated, see Figure (6). 



for two different values of the window, A. Let us stress that the difference from zero of the 
residual values here defined give a direct measurement of the level of correlation between the 
two fluctuating quantities. Positive values of the difference in (7) and (8) indicate correlation 
while negative values indicate anti-correlation. In Fig. 5 we plot the case of (no-average in time) 
A = 0, and in Fig. 6 the case of A = 9r n . As one can see, in the first case, (Fig. 5), a strong 
correlation between extreme events in the acceleration and the vorticity is observed, although 
almost no difference is detectable between the longitudinal and centripetal acceleration (right and 
left panels, respectively). On the other hand, once averaged in time, i.e. filtering out the high 
oscillation in the longitudinal acceleration, a strong pick in the acceleration- vorticity contour plot 
is visible only for the centripetal case (see left panel of Fig. 6. This is in our view a simple but 
effective way to highlight the strong correlation between the observed high tails in the acceleration 
distribution and the presence of vortex filaments in the Eulerian field. 

2 Conclusions 

We have presented results on the joint Lagrangian acceleration statistics and Eulerian vorticity 
field from DNS of fully developed turbulence. In particular we have shown the existence of strong 
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Figure 6: (Left panel). Residual value, 7£(log(a^),log(f2)), for the centripetal acceleration and 
the local enstrophy after averaging in time, over a window A = 9r rr (Right panel). The same 
but for the longitudinal acceleration, 7£(log(a; ), log(f2)). The isoline is drawn for a fixed value 
of 1Z, to guide the eyes. Notice the persistence of the strong correlation present between intense 
vorticity events and intense centripetal acceleration. On the other hand, now the longitudinal 
acceleration is not any more correlated for intense events with the enstrophy (compare the bottom 
right corners of both panels). 



correlation between intense and persistent centripetal acceleration and the presence of intense 
vortical structure. No effects similar to this one is detected for the longitudinal components of the 
acceleration. We interpret this as an evidence of particle trapping in vortex filaments. Further 
investigations are necessary in order to understand if and how this events affects the local energy 
transfer mechanism and, if it is the case, how to include the presence of these filamentary structure 
in the energy cascade phenomenology. 
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supercomputing center CINECA (Bologna, Italy) and "Centro Ricerche e Studi Enrico Fermi" 
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